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INVESTIGATION  OF  COLLISIONAL  EFFECTS  IN  THE  PLASMA 
EROSION  OPENING  SWITCH 


I .  Introduction 

Inductive  storage  devices  store  energy  in  the  form  of 
magnetic  fields  rather  than  electric  fields  and  can,  in 
principle,  be  made  considerably  smaller  and  cheaper  than  their 
capacitive  counterparts.  The  plasma  erosion  opening  switch 
(PEOS)  is  used  in  inductive  storage  applications  as  a  means  of 
pulse  compression  and  vn] tage  multiplication1 *  1 2  .  Tne  switch 
consists  of  a  low  density  plasma  ( n0  .  1012  -  1014  cm-3) 
injected  with  flow  velocity,  Vp ,  between  a  transparent  outer 
anode  of  radius  rA ,  and  a  solid  inner  cathode  of  radius  rK .  In 
experiments  on  the  Gamble  I  and  II  generators,  the  ion  species  is 
composed  mostly  of  doubly  ionized  carbon  ions,  C**.11  The  switch 
plasma  conducts  the  generator  current  for  some  period  of  time 
while  the  energy  of  the  generator  is  stored  inductively  in  the 
vacuum  transmission  line  between  the  generator  and  the  plasma. 

The  conduction  phase  ends  when  the  current  in  the  switch  reaches 
a  level  determined  by  plasma  conditions  and  the  geometry  of  the 
switch  region.12  At  this  predetermined  current,  the  plasma 
ceases  to  act  as  a  low  impedance  short  to  ground  and  current  is 
rapidly  diverted  to  the  load  (opening  phase).  This  opening  is 
accomplished  by  magnetic  insulation  of  electron  flow  in  a  gap 
formed  at  the  cathode,13  presumably  by  erosion. 

It  has  been  observed  experimentally  that  during  its 
conduction  phase,  the  PEOS  allows  magnetic  field  to  penetrate 
deeply  into  the  main  body  of  the  plasma- 13-15  Both  the  depth  of 
the  penetration  and  the  axial  width  of  the  current  channel 
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induced  in  the  plasma  are  many  times  wider  than  the  collisionless 
skin  depth.  In  fact,  the  conduction  phase  ends  and  the  switch 
opens  only  when  the  magnetic  field  has  penetrated  the  entire 
plasma.  This  phenomenon  has  provoked  a  great  deal  of  theoretical 
interest.  The  problem  is  also  of  interest  because  of  its 
application  to  the  fusion  program  and,  in  particular,  to  the  PBFA 
II  light  ion  accelerator  at  Sandia  National  Laboratory.8-16 

Earlier  analytic  work  and  simulations  of  the  PEOS  using  both 
particle  and  fluid  approaches17'22  have  shown  that,  in  the 
absence  of  anomalous  collisions,  current  and  magnetic  field 
penetrate  the  plasma  during  the  conduction  phase,  but  not  in  the 
manner  expected  from  experimental  observations.  Experimental 
data13  suggest  that:  a)  current  (and  magnetic  field)  penetrates 
axially  in  broad  channels;  b)  the  penetration  of  the  field  at  the 
anode  and  cathode  occurs  at  almost  the  same,  nearly  constant, 
rate  (the  current  flow  is  predominantly  radial);  c)  the  maximum 
switch  current  before  opening  (the  "conduction  current")  scales 
almost  linearly  both  with  the  plasma  density,  n,  and  with  the 
switch  length,  l;23  and  d)  switch  opening  occurs  at  the  cathode. 
The  above  scaling  agrees  with  that  derived  from  a  model  of  the 
switch  behavior  based  on  a  bipolar  emission  description  of 
cathode  physics  which  predicts  that  the  conduction  current  scales 
linearly  both  with  the  plasma  flux,  nVF ,  and  the  switch  area, 

2  nrK  1 . 1 2 

In  several  respects,  collisionless  simulations  differ  from 
this  experimental  picture.  First  of  all,  the  penetration  of  the 
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magnetic  field  occurs  at  different  rates  and  via  distinct 
mechanisms  near  the  cathode  and  the  anode.  Although  separate 
mechanisms  may  be  responsible  for  field  penetration  at  the  anode 
and  cathode  in  experiments,  penetration  has  always  been  observed 
to  occur  at  nearly  the  same  rate.  In  collisionless  simulations, 
ions  near  *-he  cathode  are  accelerated  in  the  axial  direction  via 
JxB  forces  and  form  a  space  charge  density  clump  moving  through 
the  plasma  and  toward  the  load  (the  "piston"  effect).19  This 
clump  of  ions  serves  to  concentrate  the  electrons  emitted  f-cm  a 
broad  area  of  the  cathode  into  one  narrow  (on  the  order  of  a  few 
skin  depths)  current  channel.  The  magnetic  field  near  the 
cathode  penetrates  more  and  more  deeply  into  the  plasma  as  the 
ion  clump  is  accelerated  toward  the  load. 

Near  the  conducting  anode,  on  the  other  hand,  the  current 
channel  migrates  toward  the  load  by  two  mechanisms,  depending  on 
proximity  to  the  anode.  In  the  body  of  the  plasma,  electrons 
forming  the  current  channel  drift  radially  by  the  E.xBg  drift  and 
toward  the  load  by  the  ErxBg  drift.  Directly  at  the  anode,  where 
the  axial  electric  field  vanishes,  electrons  can  no  longer  drift 
radially,  but  must  move  axially  toward  the  load  until  their  gyro- 
radii  become  large  enough  to  allow  them  to  cross  to  the  anode. 

In  this  manner,  magnetic  field  near  the  anode  penetrates  into  the 
plasma  by  the  axial  drift  of  the  electron  current  channel.*9 
When  the  current  channel  reaches  the  load  end  of  the  plasma,  it 
disconnects  from  the  plasma  and  the  switch  opens.  At  that  point, 
electrons  become  magnetically  insulated  and  flow  toward  the  load 
rather  than  between  the  electrodes.  In  the  regimes  considered 
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in  this  paper,  the  penetration  of  the  magnetic  field  is  faster  at 

the  anode  than  at  the  cathode,  so  the  conduction  time  is 

» 

determined  by  the  mechanism(s)  controlling  the  anode  penetration. 
Two  considerations  alter  the  anode  dominated  picture:  a)  in  a 
switch  which  is  longer  in  axial  extent  and  which  conducts  for  a 
longer  time,  it  is  possible  that  the  JxB  acceleration  of  the 
piston  at  the  cathode  could  control  the  opening  process,  and  b) 
the  anode  may  not  be  a  perfect  conductor.  In  experiments,  a  wire 
mesh  or  a  set  of  widely  spaced  conducting  rods  define  the  anode 
plane  in  the  PEOS ,  and  the  Et  field  does  not  vanish  at  all  points 
of  the  anode  plane. 

Besides  the  different  penetration  rates  at  the  electrodes, 
other  disparities  with  experimental  results  also  exist.  As 
mentioned  above,  collisionless  simulations  (both  PIC  and  fluid) 
show  consistently  narrow  current  channels  ( -  a  few  skin  depths). 
Analysis  of  this  phenomenon  and  discussion  of  some  effects  near 
the  cathode  which  may  widen  the  channels  are  given  in  Refs.  [19], 
[20],  and  [22].  In  this  paper  it  will  be  shown  that  conduction 
currents  scale  roughly  like  12/S  and  n1  / *  in  the  collisionless  PI 
simulations,  rather  than  scaling  linearly  as  observed 
experimentally.  In  addition,  the  conduction  current  in 
collisionless  simulations  will  be  seen  to  be  insensitive  to 
plasma  flux.  Since  these  simulation  results  compare  poorly  with 
experiment,  it  is  concluded  that  the  conduction  physics  in  the 
collisionless  simulations  is  different  from  the  one  operating  in 
expe  r iment . 
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It  has  long  been  speculated  that  an  anomalous  collision 
mechanism  exists  in  the  plasma  and  is  responsible  for  the  broad 
current  channel  observed  in  experiments.12'13'24'25  Beside  the 
usual  instabilities,  such  as  the  Buneman,  ion-acoustic,  etc., 
which  may  be  active  in  the  switch,  another  electrostatic 
instability  has  been  recently  identified  with  many  germane 
properties.2®  This  mode  has  a  linear  growth  time 
(Ta-7/wpt)  of  about  2-3  ns  for  typical  experimental 
parameters;  its  effective  collision  frequency  is  large  ( < 

O..4o >  )  because  it  is  weakly  saturated  and  grows  to  large 

amplitudes;  finally,  the  mode  propagates  at  the  ion  acoustic 
speed  so  that  it  can  achieve  this  large  amplitude  without 
convecting  out  of  the  plasma  during  the  time  of  interest. 
Simulations  are  currently  being  conducted  to  investigate  this 
instability  further,  especially  its  non-linear  behavior  and 
dependence  on  plasma  parameters.  To  include  this  instabi  lity 
self-consistently  into  the  current  class  of  two  (spatial) 
dimensional  particle  codes  is  impossible  because  the  unstable 
wave  requires  a  component  in  the  third  dimension.  Even  with  a 
three  dimensional  code,  a  completely  self-consistent  treatment  of 
the  problem  may  be  difficult  and  expensive.  The  dimensions 
associated  with  experimental  geometries  is  large  compared  with 
the  cell  sizes  needed  to  resolve  the  cathode  emission  process  and 
the  wavelength  of  the  unstable  modes.  A  large  number  of  macro- 
particles  per  cell  are  also  necessary  to  resolve  the  Maxwellian 
distribution  and  provide  enough  electrons  resonant  with  the 
unstable  wave.  For  these  reasons,  an  ad  hoc  method  will  be  used 
here  to  include  the  effect  of  the  instability  in  the  PIC 
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simulations  of  the  PEOS  problem.  The  use  of  a  particle  code  in 
the  r  -  0  plane  to  examine  the  instability  in  a  small  volume  will 
be  the  subject  of  a  future  paper. 

In  the  present  treatment,  a  time  varying  region  is  chosen 
where  the  instability  is  assumed  to  be  active  and  then  a 
collision  frequency,  va ,  is  selected.  The  collision  term  is 
incorporated  in  the  electron  equations  of  motion  so  that 
electrons  are  subjected  to  a  damping  force  -\>4myV,  where  y  is  the 
relativistic  factor.  The  collision  frequency  is  chosen  to  be 
within  the  range  of  values  predicted  by  the  instability  analysis, 
but  its  exact  value  and  its  region  of  influence  are  not  chosen 
self-consistently  or  with  reference  to  turn-on  conditions,  turn¬ 
off,  saturation,  etc.  The  underlying  motive  is  to  discover  what 
effect  such  an  instability  would  have  on  the  behavior  of  the 
PEOS.  Thus  the  purpose  of  this  paper  is  not  to  prove  that  such 
an  instability  exists  or  to  study  tne  exact  details  of  its 
effect,  but  to  show  that  its  presence  could  alter  the  behavior  of 
the  simulated  switches  so  that  they  operate  more  like  those 
observed  experimentally. 

The  effects  of  anomalous  resistivity  have  been  investigated 
previously  in  particle  and  fluid  simulations  of  the 
PEOS . 2 2  ■  2 5 ' 2 7 “ 2 9  Although  net  treating  the  instability  described 
above,  the  usual  streaming,  ion  acoustic,  and  lower  hybrid 
instabilities  have  been  modeled  in  detail  by  including  the 
appropriate  transport  coe f f i ci ents . 2 8 ■ 2 9  Compared  to 
collisionless  results,  these  simulations  have  shown  improved 
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agreement  with  experiment  concerning  channel  widths  and 
penetration  velocity.  In  this  paper,  a  more  general  study  of 
anomalous  collision  effects  is  made,  with  the  focus  on  scaling  of 
PEOS  operation  with  length,  density,  and  plasma  injection 
velocity  for  a  fixed  collision  rate.  Scaling  studies  are  also 
performed  for  collisionless  simulations  and  compared  with  the 
collisional  results. 


II .  Description  of  Simulations 

The  work  presented  in  this  paper  treats  the  conduction  of 
current  through  a  low  density  ( n0  -  0.25-4xl012  cm-3)  plasma 
injected  between  two  coaxial  conducting  cylinders  as  shown  in 
Fig.  1.  The  plasma  is  injected  through  the  anode  toward  the 
cathode  with  flow  velocity,  Vr ,  and  is  assumed  to  be  azimuthally 
symmetric.  The  coaxial  electrodes  and  the  switch  plasma  form 
part  of  a  circuit  with  the  generator  at  the  left  and  a  load  in 
parallel  with  the  switch  at  the  right.  In  the  present 
simulations  where  the  fCwUS  is  on  the  conduction  phase  of  the 
switch,  the  load  is  chosen  to  be  a  short  circuit.  The  regions  on 
both  sides  of  the  plasma  between  the  plasma  and  the  axial 
boundaries  are  in  vacuum.  An  electromagnetic  pulse  is  driven 
into  the  system  from  the  left  and  induces  current  in  the  plasma. 
Current  rises  to  about  120  kA  in  5  ns.  The  time  scale  of  5  ns  is 
chosen  to  reduce  the  computation  time.  To  partially  compensate 
for  this  artificially  small  time  scale,  the  ion  species  is  chosen 
to  be  hydrogen,  H*  ,  rather  than  carbon,  C**. 


Space  charge  limited  electron  emission  is  allowed  on  the 
conducting  surfaces  and  particles  that  encounter  the  surfaces  are 
absorbed.  The  2-1/2D  electromagnetic  particle-in-cell  code 
"MASK"30  is  used  to  treat  the  switch  plasma,  the  boundary 
emission  of  electrons  into  the  plasma,  and  the  initial  flow  of 
plasma  into  the  simulation  region.  The  problem  is  done  in  r-z 
geometry  and  azimuthal  symmetry  is  imposed.  Three  separate 
issues  are  addressed  here.  These  are  conduction  current  scaling 
with  length,  density,  and  plasma  flow  velocity. 

Scaling  simulations  are  run  with  and  without  anomalous 
collisions.  Anomalous  collisions  are  introduced  in  regions  of 
the  plasma  in  which  the  magnetic  field  is  above  a  certain 
threshold  (  10%  of  the  maximum  B  field  at  any  given  time)  by 
turning  on  the  collisional  term  -v4 rnyV  in  the  particle  equations 
of  motion.  The  magnetic  field  is  used  as  a  turn  on  condition 
because  its  presence  signals  that  current  is  being  driven  through 
the  plasma  which,  in  turn,  can  drive  the  instability.  As 
mentioned  earlier,  the  motivation  here  is  to  determine  whether 
collisional  simulations  of  the  switch  more  closely  resemble 
experimentally  observed  scalings. 

Although  the  length,  density,  and  plasma  flow  velocity  are 
varied  in  the  simulations  discussed  herein,  the  same  driving 
current  waveform  (rising  linearly  to  about  120  kA  in  5  ns;,  the 
same  radial  dimensions  ( rK  =2.5  cm,  =  5  cm),  and  the  same 
collision  frequency  (in  the  collisional  runs),  \>a  =  i.GxlO15  s' :  , 
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ace  used  in  all  cases.  There  was  no  attempt  to  scale  \>,  with 
density,  magnetic  field,  or  other  parameter  as  dictated  by  any 
given  instability.  However,  for  reference  purposes,  note  that 
the  ratio  varies  from  0.09  for  no  ■  4x10* 2  cm* 3  to  0.35 

for  n0  -  2.5X1011  cm- 3  • 


III.  Results 


This  section  begins  with  an  examination  of  length  and 
density  scaling  for  the  series  of  simulations  described  in  Table 
I.  The  discussion  is  broken  into  two  sections,  one  devoted  to 
collisionless,  the  other  to  collisional  results.  These  two 
sections  are  followed  by  a  third  concerning  the  effects  of 
varying  the  injection  velocity, 

A.  Collisionless  Simulations 


In  Fig.  2,  the  motion  of  the  current  front,  D(t),  is  plotted 
as  a  function  of  time  for  the  case  where  n  -  1C12  cm* 3  ,  1  *  10  cm 
and  Vr  -  1.8x10*  cm/s.  The  position  of  the  current  front  near 
the  anode  (r  -  4.5  cm)  and  near  the  cathode  (r  »  3.0  cm)  are 
plotted  separately.  The  current  front  position  is  defined  to  be 
the  axial  position  where  B 9  begins  to  rise  from  zero.  A  least 
squares  power  curve  fit  is  used  on  the  data  points.  The  current 
penetration  distance  during  the  first  0.4  ns  is  not  represented 
accurately  in  these  curves  because  the  initial  penetration,  5, 
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occurs  rapidly  over  the  distance  of  a  few  skin  depths,  6  =  2.4 
cm. 

The  anode  and  cathode  curve  fits,  Dx  ( t )  and  DK ( t ) ,  are 

tabulated  in  Table  I,  case  (a).  The  anode  front  scales  as  t2  ■  4 , 

while  the  cathode  front  as  t2 • 5  .  Since  the  current  rises 
linearly  in  time,  and  the  switch  opens  in  this  simulation  when 
the  anode  front  reaches  the  load  end  of  the  plasma  <Dx(t)  -  1), 
it  is  inferred  that  the  conduction  current  scales  with  length 
roughly  as  11/2-*.  Similar  scaling  can  be  done  for  all  the  other 
collisionless  runs  shown  in  Table  I  (cases  (b)  -  (e)).  Taking  a 
rough  average  of  all  the  exponents  of  time  in  the  function  Dx  ( t )  , 
the  conduction  current  scales  roughly  as  12/5  with  length  in  the 
collisionless  case. 

In  Fig.  3,  the  current  front  position  is  shown  for  the  case 

of  a  higher  density  plasma,  n  ■  4xl012  cm-3  ,1-10  cm,  and  vF  = 

1.8x10®  cm/s .  The  initial  penetration  is  reduced  to  6  =  1.1  cm, 
and  the  anode  penetration  curve  scales  as  t2 • 3  with  time,  as  seen 
in  Table  I,  case  (b).  The  functional  dependence  of  the  cathode 
end  of  the  current  channel  on  time  varies  strongly  as  the  density- 
changes  from  1012  cm- 3  to  4xl02 2  cm- 3  .  As  shown  in  Table  I, 

DK ( t )  -  0.5t2-5  for  n0  -  1012  cm"3,  while  DK ( t )  -  O.lSt1-4  for 
n0  -  4xl012  cm" 3  .  Cathode  penetration  for  the  n0  =  4xl012  cm" 3 
case  follows  the  axial  penetration  of  a  leaky  piston  accelerated 
by  JxB  forces  through  the  plasma.19  However,  at  nc  = 

1012  cm"3,  the  current  front  position  and  piston  do  not  track 
because  a  small  amount  of  current  is  able  to  penetrate  4-5  cm 


ahead  (toward  the  load)  of  the  piston.  The  amount  of  current 
more  than  1  cm  in  front  of  the  piston  is  less  than  30%  of  the 
total  carried  by  the  switch,  and  the  amount  of  current  more  than 
2  cm  in  front  is  less  than  15%  of  the  total.  So  the  piston  still 
defines  the  most  active  region  of  electron  emission.  The  piston 
in  the  n0  -  1012  cm- 3  run  depends  on  time  roughly  as  the  function 
0 . 4  tl  •  4  . 

Several  other  regimes  of  density  and  length  have  been 
examined  and  results  are  shown  in  Table  I,  together  with 
estimates  of  the  scaling  of  conduction  current  with  density.  All 
cases  in  Table  I  have  Vr  -  1.8x10*  cm/s.  Comparing  the  two  1-5 
cm  cases  ((d)  and  (e)),  at  nQ  -  1012  cm- 3  and  4xl012  cm-3,  we 
find  from  the  actual  time  when. opening  begins,  x0 ,  that  the 
conduction  current  scales  as  n1/*.  The  density  scaling  is 
obtained  by  assuming  that  I0  -  x0  ~  na  for  a  linearly  rising 
current  drive.  When  a  correction  is  made  for  the  different 
initial  penetration  depths  (2.4  cm  and  1.1  cm,  respectively),  the 
conduction  time  scales  more  like  nl/1-  The  corrected  scaling  is 
obtained  by  assuming  the  same  initial  penetration  depth  (2.4  cm 
for  (a),  (b)  and  (c),  2.5  cm  for  (d)  and  (e))  applies  to  the 
cases  to  be  compared.  A  new  estimated  opening  time  is  found  by 
solving  Da(t)  ■  1  for  t0  (with  the  modified  6),  and  then  used  to 
scale.  This  scaling  more  accurately  reflects  the  actual 
penetration  velocity  through  the  plasma  after  the  initial  current 
penetration  has  occured.  Corrected  scalings  for  1  =  10  cm 
comparing  n0  =  101 2 ,  4xl012,  and  2.5xl013  cm- 3  show  scaling  as 


The  conclusions  reached  above  for  length  and  density  scaling 
pertain  to  collisionless  switch  simulations  for  which  anode 
current  penetration  dominates  cathode  penetration.  Collisionless 
simulations  of  the  switch  using  the  implicit  code  "ANTHEM"  in  the 
parameter  regimes  of  Gambles  I  and  II  (x  =  50  ns,  C+*  plasma) 
with  a  conducting  anode  also  show  opening  dominated  by  anode 
penetration,29  so  the  conclusions  obtained  here  for  a  t  a  3-5  ns, 
H*  plasma  should  apply  to  the  Gamble  regimes.  For  a  longer 
conduction  time  switch  (i.e.,  x0  >  100  ns  with  a  C* *  carbon 
plasma),  JxB  acceleration  of  the  ion  piston  could  dominate  the 
opening  process  in  a  collisionless  simulation  of  the  switch  and 
stronger  scaling  with  density  and  length  would  be  seen.  The  same 
remark  applies  to  the  case  when  the  anode  is  not  a  perfect 
conductor.  If  cathode  penetration  ( 0K ( t ) )  determined  the  opening 
time  in  these  simulations,  one  would  find  nearly  linear  corrected 
scaling  with  density  when  comparing  n0  -  1012  cm* 3  and  n0  * 

4xl012  cm* 3  (with  6  -  2.4  cm,  1-10  cm).  However,  when 
comparing  n0  -  2.5xlOxl  cm* 3  and  n0  -  1012  cm* 3  ,  the  scaling 
becomes  n° • 4 .  This  change  in  the  scaling  may  be  the  result  of 
the  effect  described  earlier,  namely  that  the  cathode  current 
front  tracks  the  piston  at  high  densities  but  can  move  slightly 
ahead  of  the  piston  at  low  densities.  Conduction  current  scaling 
with  length  can  also  be  inferred  from  the  indices  of  time  in  the 
functions  DK ( t ) .  Again  the  scaling  for  the  1  =  10  cm  runs  seems 
to  be  a  function  of  the  density  regime.  At  densities  of  10;; 
cm* 3  or  lower,  the  conduction  current  would  scale  roughly  as 
11/2-3,  while  at  n0  *  4xlOx 2  cm* 3  it  scales  as  l171-9. 
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Summarizing  the  collisionless  results,  it  is  concluded  that 
conduction  current  scales  roughly  as  12/5  with  length  and  as  n1/4 
with  density  in  this  anode  dominated  regime.  Stronger  scalings 
are  found  from  the  cathode  penetration  distances.  However,  the 
scaling  is  dependent  on  the  density  regime.  At  the  cathode, 
current  penetrates  primarily  by  the  axial  JxB  acceleration  of  a 
leaky  ion  piston,  while  at  the  anode  axial  ExB  drift  of  electrons 
moves  the  front  toward  the  load.  The  current  channels  observed 
in  collisionless  runs  are  usually  about  1  cm  wide  for  no  -  4xl0: 2 
cm-3,  1.5-2. 5  cm  wide  for  n0  -  1012  cm"3,  and  3-4.5  cm  for  n0  - 
2.5xl013  cm-3.  These  widths  are  about  3-4  skin  depths  at  the 
various  densities. 

The  anode  penetration  distances  presented  here  depend  on 
time  to  some  power  greater  than  two,  while  similar  distances 
presented  in  Ref.  [19]  depend  logarithmically.  This  apparent 
discrepancy  is  a  result  of  the  radial  location  of  the 
measurement.  In  the  present  paper,  the  axial  penetration 
distance  is  measured  at  r  -  4.5  cm.  In  the  previous  work  it  was 
taken  at  r  »  5.0  cm  where  the  axial  ExB  drift  dominates  because 
Et  and  the  radial  ExB  drift  vanish  at  the  conducting  anode. 

These  differences  again  give  an  indication  of  how  important  the 
anode  boundary  conditions  are. 

B .  Collisional  Simulations 

Current  penetration  distances  from  three  collisional 
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simulations  are  shown  in  Fig.  4.  The  parameters  used  are 
outlined  in  Table  I,  cases  (f)-(h).  Only  the  first  3.2  ns  of  the 
n  •  4xl012  cm" 3  run  are  shown.  In  the  lower  density  regimes. 

O 

the  anode  and  cathode  currents  migrate  through  the  plasma  at 
approximately  the  same  rate,  so  only  one  curve  is  drawn.  At 
higher  densities  (>  2xl03 2  cm-3),  the  penetration  rates  are 
different.  This  occurs  because  the  collisional  frequency  v4  - 
1010  s" 1  used  in  the  higher  density  cases  is  not  sufficient  to 
prevent  the  JxB  piston  effect  described  in  the  introduction. 

When  a  higher  collision  frequency,  v4  -  2xl010  s_1,  is  used  in 
the  ne  -  4xl01J  cm" 3  case,  anode  penetration  remains  almost  the 
same  as  for  v4  -  1010  s" 1 ,  but  cathode  penetration  proceeds  at  a 
faster  pace.  At  t  -  3.2  ns  with  \>4  -  2xl010  s“ 1 ,  the  cathode 
current  channel  is  0.8  cm  behind  the  anode  current  channel,  while 
with  v4  -  1010  s" 1 ,  the  channels  are  separated  axially  by  1.3  cm. 
That  collisions  have  a  more  substantial  effect  at  the  cathode 
than  at  the  anode  is  consistent  with  earlier  conjecture  that 
anomalous  collisions  at  the  cathode  alone  may  be  sufficient  to 
cause  broad  channels  and  a  deeply  penetrating  magnetic 
field .22.25 

The  first  observation  from  Fig.  4  is  that  the  penetration  is 
linear  in  time  rather  than  a  power  of  time.  The  conduction 
current  therefore  scales  linearly  with  switch  length.  Figure  4 
also  shows  that  the  rate  of  penetration  is  more  strongly 
dependent  on  the  density  than  in  the  coilisionless  simulations. 

As  detailed  in  Table  I,  the  slope  of  the  anode  curve  is  2.9  c m/ns 
at  n0  -  1012  cm-3,  while  it  is  1.34  cm/ns  at  n0  =  4xlOi;  cm" 3 
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and  2.2  cm/ns  at  ne  -  2xlOiJ  cm- 5 .  At  the  lower  density,  n0  - 
2.5X1011  cm-3,  the  slope  is  4.9  cm/ns.  The  scalings  for  these 
cases  show  that  the  conduction  time  scales  roughly  as  n1/2  in  the 
anode  dominated  case. 

If  the  conduction  process  were  cathode  dominated  (because 
the  anode  is  not  a  perfect  conductor,  for  example),  then  even 
stronger  scalings  would  apply.  According  to  the  DK ( t )  curves  in 
the  collisional  case,  conduction  current  scaling  is  linear  in 
both  density  and  length.  In  contrast  to  the  collisionless  case, 
these  scalings  are  independent  of  density  regime.  The  mechanism 
behind  cathode  penetration  in  the  collisional  simulations  is  no 
longer  axial  JxB  acceleration  as  in  the  collisionless  cases.  The 
scaling  of  cathode  penetration  with  density  and  length  suggests 
that  the  mechanism  is  related  to  bipolar  space  charge  limited 
emission  as  described  in.12 

An  additional  important  observation  concerning  the 
collisional  results  is  that,  as  found  in  previous  work,  the 
current  channels  are  wider  than  in  the  collisionless  regime.  It 
is  found  that  for  a  10  cm  plasma,  channel  widths  of  7-9  cm  are 
typical.  The  current  is  spread  fairly  uniformly  early  in  a 
simulation,  but  becomes  somewhat  more  concentrated  near  the 
generator  side  of  the  plasma  late  in  time  after  plasma  density 
accumulates  in  a  weak  ion  piston  which  forms  near  the  cathode. 

In  Figs.  5  -  7,  a  comparison  is  made  between  collisionless  and 
collisional  results.  The  figures  are  taken  from  cases  (a)  and 
(f)  of  Table  I.  The  axial  length  of  the  region  shown  is  12.8  cm, 
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and  its  radial  height  is  2.5  cm.  Dashed  vertical  lines  in  the 
simulation  region  define  the  boundary  of  the  original  plasma 
fill.  Figure  5  shows  magnetic  field  contours  at  t  -  2.4  ns  for 
the  collisionless  and  the  collisional  case.  The  contour  levels 
are  the  same  in  both  plots  and  vary  from  -4.5  kG  (level  1)  to  - 
55.56  G  (level  9)  with  an  increment  of  0.55  kG  per  level.  The 
figure  shows  that  the  magnetic  field  is  spread  uniformly  over  a 
width  of  about  7.6  cm  in  the  collisional  simulations,  with  a 
slight  concentration  of  current  near  the  generator  side  of  the 
plasma.  In  the  collisionless  simulations,  current  is  carried 
diagonally  between  the  electrodes  in  a  channel  of  width  -  2  cm. 

In  Figs.  6  and  7,  electron  and  ion  particle  plots  compare 
results  from  the  collisionless  and  collisional  simulations  at  t  * 
2.4  ns.  As  reported  in  earlier  work,19  in  the  collisionless  case 
electrons  E::3  drift  diagonally  between  the  electrodes,  with  the 
electric  field  in  the  drift  channel  produced  largely  by  net 
positive  space  charge  to  the  left  of  the  diagonal  electron 
current  channel.  Ion  charge  accumulates  near  the  cathode  in  a 
diagonal  piston  which  JxB  accelerates  axially  toward  the  load. 

In  contrast,  the  electron  flow  in  the  collisional  case  (see 
Figs.  5  (b)  and  6(b)),  is  predominantly  radial  and  spread 
throughout  the  plasma.  Certain  effects  from  the  collisionless 
regime  are  also  manifest.  For  example,  diagonal  flow  is 
established  at  the  generator  edge  of  the  plasma,  but  some 
electrons  spiral  in  an  ExB  drift  which  does  not  allow  them  to 
cross  to  the  anode.  That  no  current  crosses  to  the  anode  at  the 
generator  side  of  the  plasma  can  be  seen  by  studying  Figs.  5(b) 
and  6(b).  As  in  the  collisionless  regime,  the  electric  field 
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responsible  for  the  ExB  drift  is  produced  by  space  charge 
separation.  Also,  a  piston  forms  in  the  ions  even  in  the 
collisional  case,  although  generally  it  is  smaller  and  takes  a 
few  nanoseconds  longer  to  appear. 


Plasma  Iniection  Velocity  Effects 


An  examination  has  been  made  of  the  sensitivity  of  the 
current  penetration  velocity  on  the  plasma  injection  velocity, 

Vr .  Although  no  direct  experimental  evidence  exists  for  this 
scaling,  the  NRL  theory12  of  current  penetration  velocity  scaling 
with  vf  suggests  that  it  should  be  linear,  with  a  correction 
factor  depending  on  the  ratio  of  the  cathode  gap  potential,  e$, 
to  the  ion  injection  energy  MiVf2/2.20  The  space  charge  limited 
current  density  is  given  by, 


'bp  ' 


where 


2Ze  ♦ 

M.  V  ‘ 
i  F 


Here  Z  is  the  ion  charge  state  and  Mt  the  ion  mass.  In  Eq.(l), 
jbp  is  the  usual  bipolar  current  density  relating  the  electron 
current  density  to  the  smaller  ion  current  density,  n,  Ievr ,  so 
that , 
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Note  that  the  correction  factor,  a,  results  from  injecting  ions 
into  the  accelerating  gap  with  non-zero  energy.  From  Eqs .  (D- 
(4),  a  1012  cm-3,  H*  plasma  of  length  1  ■  5  cm,  injected  at  Vr  - 
1.8x10*  cm/s  into  the  system  depicted  in  Fig.  1  is  expected  to 
conduct  approximately  I0  -  65  kA  before  opening  (assuming  that  nt 
-  1012  cm" 3  at  the  cathode  and  that  the  gap  potential  *  «  10s  v, 
as  is  typical  in  the  simulations).  In  the  collisionless  and 
collisional  simulations,  however,  I0  =  37  kA,  suggesting  that 
other  factors  also  affect  the  conduction  current  I0  .  Scaling 
with  Vr  can  still  be  tested  by  assuming  that  IQ  =  37  kA  is  the 
space  charge  limited  current  for  Vr 0  -  1.8x10*  cm/s  and  then 
injecting  Vr  -  4Vro  and  Vr  .  vro/4  to  see  Io  increases  or 
decreases  incrementally  as  expected  from  Eqs.  (l)-(4).  The 
results  are  shown  in  Table  II.  From  left  to  right  in  the  table 
are  the  collision  frequency,  injection  velocity,  cathode  gap 
potential,  the  factor  a  from  Eq .  (2),  and  the  predicted 
conduction  current,  I,p.  The  current  Isp  is  found  by  using  Eqs. 

( 1 ) — { 4 )  and  then  adding  the  ion  current  at  the  cathode  in  order 
to  obtain  the  full  conduction  current,  Isp  =  Ie  +  lL.  Because  rK 
■  rA/2  and  because  the  initial  plasma  density  is  uniform  in  r, 
the  cathode  ion  current  is  one  half  the  ion  current  injected  at 
the  anode,  IF ,  shown  in  the  next  column  in  Table  II.  A 
comparison  of  I5p  to  the  actual  conduction  current,  I0  ,  in  the 
right-most  column  shows  poor  cor respondence  in  both  the 
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collisional  and  collisionless  regimes.  Little  change  is 
observed  in  Ia  as  Vr  is  varied  in  the  collisionless  regime. 
Collisional  results  do  show  sensitivity  to  Vf ,  but  the  scaling  is 
weaker  than  expected  from  Eqs.  ( 1 ) — ( 4 ) . 

In  the  collisional  runs,  a  curious  correspondence  exists 
between  the  additional  conduction  current  carried  by  the 
simulated  PEOS  and  the  additional  ion  current,  Ip  ,  injected  at 
the  anode  into  the  system.  Of  course,  both  electrons  and  ions 
are  injectad  at  the  anode  with  the  same  injection  velocity,  but 
many  electrons  begin  to  ExB  drift  axially  so  that  a  radial  ion 
current  Ip  -  2nlrAn1ZeVp  becomes  available  near  the  anode.  A 
simulation  of  a  stationary  plasma  (Vr  -  0)  shows  that  it  carries 
conduction  current  I,  -  33  kA.  Using  this  as  a  base  current,  the 
conduction  current  due  to  the  additional  ion  current,  lr ,  becomes 

shown  in  the  next  to  last  column  (1^  -  If  +  Ir  )  .  In  the 
collisional  regime,  IAP  predicts  the  conduction  current  fairly 
well  for  all  injection  velocities  except  the  1.44x10*  cm/s  case. 
It  is  possible  that  the  (initially  comoving)  electron  current  is 
not  deflected  axially  as  effectively  at  these  high  velocities. 

The  conclusions  reached  above  apply  to  the  case  simulated 
above  where  the  anode  is  a  perfect  conductor.  If  the  anode 
allowed  an  axial  electric  field,  electrons  could  drift  radially 
near  the  anode  and  different  scaling  with  VF  may  be  seen. 


IV. 


Conclusion 


An  attempt  has  been  made  in  this  paper  to  reconcile 
experimental  observations  and  numerical  results  regarding  the 
conduction  phase  of  the  PEOS .  Experimentally  it  is  found  that 
current  and  field  penetrate  axially  into  the  plasma  in  wide 
channels.13'15  More  recently,  the  conduction  current  has  been 
observed  to  scale  roughly  linearly  with  plasma  length  and 
density.23  An  analytic  model  of  the  switch  predicts  that  the 
conduction  current  should  scale  linearly  with  density,  length, 
and  injection  velocity.12  In  contrast  to  these  expectations, 
collisionless  simulations  of  the  PEOS  in  the  anode  dominated  case 
show  narrow  channels,  r\l  /  *  dependence  on  density,  almost  no 
sensitivity  to  injection  velocity,  and  12/5  scaling  of  the 
conduction  current  with  length.  From  an  examination  of  the 
cathode  penetration  curves  in  the  collisionless  simulations, 
conduction  current  scalings  with  length  and  density  depend  on  the 
density  regime  in  question.  Generally,  the  scalings  are  stronger 
at  the  cathode  than  the  anode,  but  still  sublinear.  The  current 
penetrates  at  the  cathode  by  convective  JxB  acceleration  of  a 
piston  toward  the  load,  although  at  low  densities  a  small  amount 
of  current  is  able  to  penetrate  ahead  of  the  piston. 

When  an  ad  hoc  collisional  model  is  introduced  in  the 
simulations,  agreement  with  experiment  improves.  In  these  anode 
dominated  simulations,  conduction  current  scales  linearly  with 
length,  as  n1/2  with  density,  and  much  wider  current  channels  are 
produced.  Although  some  sensitivity  to  injection  velocity  is 
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seen,  the  dependence  is  much  weaker  than  predicted  from  a  space 
charge  limited  model  of  the  conduction  phase.  Even  better 
agreement  with  experiment  is  obtained  from  the  cathode 
penetration  curves  in  the  collisional  simulations.  It  can  be 
inferred  that  that  if  the  conduction  phase  of  the  switch  were 
cathode  dominated,  then  the  conduction  current  would  scale 
linearly  with  both  length  and  density.  Unlike  the  collisionless 
results,  these  scalings  are  independent  of  the  density  regime 
being  considered.  The  cathode  penetration  mechanism  in  the 
collisional  case  is  no  longer  axial  JxB  accele ra t i on ,  and  seems 
consistent  with  the  NRL  model.12 

Future  work  in  this  area  should  confirm  the  cathode 
dominated  scalings  using  non-conducting  anode  boundary 
conditions.  More  work  is  also  needed  to  examine  the  effects  of 
improved  modeling  of  instabilities,  and  to  understand  the  weak 
conduction  current  scaling  with  injection  velocity. 
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Table  I.  Anode  ( Dx )  and  cathode  ( DK )  penetration  distances  and 
conduction  times  (t. )  for  various  collisionless  and  collisional 


runs . 


t/  X  1010  n  x  1012  1 

DA(t ) 

DK(t) 

(j“l)  (cm'3)  (cm) 

(cm) 

(cm) 

2.4  +  0.66124 

1.1  +  0.34123 

4.4  -I-  2.412-2 

2.5  +  0.64f2  4 

1.2  +  0.31<2  7 

1.3  +  2.9f 
0.7  +  1.341 

3.3  +  4 .91 

1.1  +  2.2f 


2.4  -f  0.47l2s 

1.1  +  0.26f 1 4 

4.4  4-  2.2 12  2 

2.5  +  0.5612 

1.2  +  0.2 1 2 1  8 

1.3  +  2.9 1 

1.3  +  0.741 

3.3  +  4.9' 
1.37  +  1.51 
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Table  II.  Conduction  currents  ( I0 )  at  different  injection 
velocities  (Vr )  for  various  collisionless  and  collisional  runs 


run 

x  1010 
(*“*) 

a 

Effl 

h 

(IcA) 

If 

(kA) 

a 

0 

1.8  1  85 

0.65 

- 

KOI 

-  |  37.5 

mm 

0 

7.2  |  89 

0.27 

64.5 

18.1 

51.1 

39.0 

c 

1.0 

1.8 

100 

0.67 

- 

MEM 

Kfl 

1.0 

7.2 

100 

0.28 

66.7 

18.1 

51.1 

50 

n 

1.0 

0.45 

77 

0.89 

12.5 

1.13 

ItEMI 

f 

1.0 

14.4 

100 

0.15 

76.2 

36.2 
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Fig.  1.  Simulation  region  showing  components  of  typical  PEOS 
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p)  and  a  coliisional  (run(f)  at  bottom)  case  at  t  =  2.4  ns. 
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